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ABSTRACT 


This work presents comprehensive study of MOSFET model. Also presented are models 
lor polysilicon TFT in subthreshold and accumulation regions. The MOSFET model studied 
enables greater understanding of relatively short-channel (> 1 fj,m) devices in terms of drain 
current vs drain voltage characteristics. This model takes into account effect of surface scat- 
tering and channel scattering on the carrier mobility. The parameters 0 O and 6\ arc used to 
model effects of surface scattering and intrinsic resistance of source and drain region on the 
carrier mobility, the* effects of the parameters % and 771 account for the effects of lateral field 
in the channel on the carrier mobility. These effects include the channel velocity saturation. 
Also present ed is a discussion on the extraction of these parameters. The influence of bulk 
doping on the drain current is described in terms of parameter ’a’. 

The subthreshold model for TFT assumes exponential density of states for the gap and ex- 
pression for differential channel conductance is described and G-Vg characteristics obtained. 
Tlu* model for accumulation region is based on the assumption that the potential barrier 
height at grain boundary depends on the gate and drain voltage of polysilicon TFT. An 
effective channel mobility with contribution from grain mobility and boundary mobility is 
, defined and In ~ 1 7> characteristics of polysilicon TFT is presented. 
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Layout Of The Thesis 


Chapter. 1. presents introduction of the thesis. 

chapter 2 describes theory of MOSFET model including derivations of model equ 
In chapter 3. there is discussion of accurate yet simple MOSFET model for enhai 
mode and depletion mode devices. This chapter uses equations which have been de 
chapter 2. Chapter 2 also presents theoretical study of MOSFET model. 

Model for subthreshold characteristics of polysilicon Thin-Film Transistors has b< 
sented in chapter 4. 

In chapter 5 modelling of polysilicon Thin-Film-Transistors in accumulation rt 
presented. 

Appendix B presents general mathematical formulations for poly-Si TFT. 



Chapter 1 




Introduction 


Over the years modelling of MOSFET has been used to facilitate the understanding o: 
MOS device. Modelling also enables to consider the various complex effects ,e.g. short c 
nel effects, influencing the performance cf the device. In the analytical models considers 
of these effects often becomes unwieldy and assumptions are made to make equations sin 
such as the assumptions that field is much larger than the lateral field. These assump 
result in low degree of accuracy in predicting device the characteristics. 

In this work, I have considered MOSFET dc model including short channel effects, 
ther, I have also considered modelling of polysilicon thin film transistor (TFT) whit 
emerging as an important new technology for large area microelectronics in liquid cr 
displays and imaging arrays. 

1.1 Scope of the Present Work 

In this work, modelling of crystalline silicon MOSFET and polycrystalline silicon TFT 
been considered. The commonly used analytical model used by circuit designers for pre 
ing MOSFET performance is based on the gradual channel approximation. 

This analytical model along with underlying assumption is presented in chapter ! 
MOSFET dc model in which several of these assumptions are eliminated is presented 
short channel device resulting in improvement in accuracy in MOSFET modelling. ; 
presented MOSFET model suitable for implementation in circuit analysis has emphasi 
carrier mobility expresssion which includes the effect of surface scattering, channel scatti 
and substrate bias which include the effects of intrinsic source and drain series resistance 
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CHAPTER 1. INTRODUCTION 


carrier mobility are also included. Polysilicon Thin Film Transistors (TFT) are emergi 
important new technology for large area microelectronics due to the higher achievable c. 
mobility in polvsilicon in comparison to that in a-Si:H. In this work modelling of subthre 
characteristics of polysilicon TFT is also presented along with model for accumulation re 
Polysicon TFT based readout electronics has potential application in solid state radi 
detector for high energy physics and for medical imaging applications. 

1.2 Statement of the purpose 

Extensive research works have been done in the area of developing analytical formula, bu 
developed formula is not capable to describe accurate and efficient short channel MOS 
modelling for Ip-V# characteristics in the entire region, subthreshold characteristics of 
Si TFT and accumulation mode of poly-Si TFT . In this thesis, a MOSFET dc n 
and TFT model is described. In which several modifications have been done to prov 
significant improvement in accuracy in MOSFET modelling for short channel and pc 
TFT modelling. The presented MOSFET model and poly-Si TFT model is very simple \ 
makes it suitable for implementation in circuit analysis programs. There is discussion 
MOSFET model with emphasis on a carrier mobility expression which includes the e 
of surface scat tering. channel scattering, and substrate bias. The effects of intrinsic si 
and drain series resistance on the carrier mobility is also included. 



Chapter 2 


Theory of MOSFET Model 


2.1 Current- Volt age Relationship 0 


Let. us consider a surface field effect transistor operating in the linear region. The voltag 
drop across an elemental section of the channel is given by, 


dV = l D dR = 


l pdy 

If fJ- oQn{y ) 


( 2.1 


Where W is the channel width and Q n (y) charge density in channel per unit area. At 


L >- 



-j- • 

The elemental section of the channel employed In the derivation of the 
current-voltage characteristics of surface field-effect transistors. 


Figure 2.1° 

distance y from the source, the total charge induced in the silicon Qs will partly consist o 
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charge in the inversion layer Q„ and partly of charge in the surface depletion region due 
the ionized acceptor ions, Qb ■ Thus, 

Qs(y) = Qn{y) + Qb(v)- (2.' 


and, 


Vg ~ V F u = - + <Ps 

(-'0 


( 2 .; 


which relates the voltage applied :o the field plate of an MOS structure to the charge induce 
in the silicon Q s and to the surface potential <f> s which denotes the total bending of the energ 
bands. We also include here the flat band voltage V FB which is due partly to the presence < 
charges in the insulator or at the interface, and partly to a finite metal-semiconductor wor 
function difference as discussed in eq.(2.3). C 0 — is the capacitance of the oxide laye 
per unit area. 

Combining these equations (2.2) & (2.3), yields the charge density in the inversion layer, 


Qn(y) = -[V G - V F b - 0,(i/)]C. ~ Qb(v) (2-4 


since it. is assumed that a conducting inversion layer exists, the surface potential 0 S wii 
be given approximately by the condition of strong inversion in the presence of an applici 
voltage. 

My) = V(y) + 2<t> f (2.5 

Whore 0/ is the Fermi potential of the substrate, V(y) is the reverse bias between t,h 

elemental section of the channel and the substrate. The charge within the surface depletioi 

! 

region Q a is given by 

Qh(v) = -Q^A-Limaxiv) = -\/(2 e s e 0 qN A [V(y) + 20/]), (2.6 1 


Since the voltage V(y) increases from the source towards the drain due to the IR drop alonj 
the channel, the field-induced junction between the n-type inversion layer and the substrata 
becomes incn'asingly r(!verse biased as it will proceed from source to drain. Hence, both tin 
energy band bonding </), and the charge within the surface depletion region Q n increase fron* 
source. to drain. Combining above equations 2.1, 2.4, 2.5, 2.6 and integrating equation (2.1^ 
between the source', where y = 0 and V=0, and the drain, where v=L and!' = Vp, yields : 


W 


< n 


I'nCoilVa - Vpn 


20 / " -^\Vd 


2 y[{2t t toqN A ) 


Co 


[(Vb + 20,) 3/2 -(20/) 3/2 ]} (2.7) 
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Which can he written in the form 


ID 


ir 

T 


v, 


ToC 0 {[V c — Vfb — 20/ — — ] Vp 




Where 


r\ J ' ls Q 

" d€ox%dmcLX,v 

~ V { Wa 


\/(20/) [{Vo + 20/) 3/2 — (20/ ) 3// “] } (2. 


) 


( 2 .: 


1 his equation (2.0) is valid only below saturation. Thus the curves shown were calculate 


(or 0 < V/> > V /w Beyond V /Xs „, the current was taken to be constant. 

It is interesting to consider the two limiting forms of Equation (2.8). For very small drai 
voltages, i.e., V,j « 20/, an expansion of the bracketed terms leads to the formulas, 


Id * ~fi 0 C 0 [ V G - V FB - 20/ + %*]Vb (2. 1C 

L 0 0 

Where 

Qb,o - -\[{2e s e 0 qN A (2(p f )) ( 2.11 

is the charge density per unit area within the surface depletion region, in equilibrium. Tliu; 
an ohmic characteristic results. The channel conductance 
in the linear region is 


fl = ~fi„C n (V(; — V r ) [linear region ] (2.12) 

Where the threshold voltage Vr, i.e., the voltage that must be applied to the gate in ordei 
to induce a conducting channel, is given by 

V r = Vfb + 20/ — —— (2.13) 

^ O 

The channel conductance of the MOS transistor in the linear region is a function of gate 
voltage. The reduction in slope .„t high gate voltages is due to a decrease in the mobility 
of electrons owing to increased surface scattering. The above derivation loses its validity 
when the inversion layer disappears near the drain, i.e., When V D — V Dsat . The condition 
Q„[L) = 0. yields 


b/w + 


s C. 0 (jNa [ViXial + 20/]) 




(2.14) 


Where if have also used equations 2.5' and 2.6, evaluated at V - V Dsat . Solving Equation 
for \ leads to 


= be; - 1/ 


i'll 


0 , , 


ci 


T 


1 + 


2 C 2 0 (Vc 


Vfb) 


c s e n qN A 


)] 


( 2 . 15 ) 
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or 


V Dsat = V G - V F3 - 2 4>f + ) 2 [1 - v/(l + (2 

€ox%dmax,o " 0 / 

Note that when the oxide thickness rr 0 is small in comparison to the width of the sur. 
depletion region x dmaXi0 , this expression reduces to the simple form, 

kjDsai l / 7 B 2 (t)f \Xo 'Edmax,o\ • (2 

subsitution of V Dsat into the current voltage relationship, This equation, gives the magnit 
of the saturation current, Iosat ■ 
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2.2 Model Equations 

Equation (2.7) re-written and modified for body bias has the form 

Id = ? ^C 0 {[V G - V FB - 2 4>, - *?\V B - + _ v 4 ,)W - - H,) 3 ' 2 ]} 

Where 7 0 = n 0 C 0 = 

and A' = . 

In the derivation of the models of MOS transistors, following assumptoins are generally made 

1. The gradual channel approximation applies. 

2. The carrier mobility in the inversion layer is constant. 

3. The reverse-biased junction leakage current is negligibly small. 

4. The drain current consists of drift current, diffusion current is ignored. 

5. In the linear region, the drain-to-source voltage is small compared to 2 <pf. 

6. Beyond saturation , the drain current is constant. 

These assumptions result in a high degree of inaccuracy in predicting current characteristics 
in short channel devices. It has been found that the following modifications of the above 
assumptions provide a significant improvement in modelling of these devices. 

1) Retention of the entire bulk doping term. 

2) Inclusion of the dependence of mobility on gate and drain voltages. 

3) Inclusion of intrinsic source and drain series resistance. 

4) Inclusion of channel length modulation. 

These four effects will be considered in the following sections. 


2.3 Bulk Doping Term 

The following equation (2.18) for linear region current results from the assumptions previ- 
ously stated, 

b = “ Vfb ~ 2 ~ ~ + 20/ " Vbs)V2 ~ (2 ^ ~ Hs)3/2] (2 ' 18) i 

It is clear that above equation reduces to 

l d = — Vt — ^2* ]^ds - i 

for V* << 2 <f>f - V bs For V ds « 2<p f - V bs , the simple model extends this approximation i 
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losult.s iii an improved fit to the measured curves, particularly as the saturation region is 
approached. To retain the bulk doping term, while including the device threshold voltage in 
an explicit form lor all values of V rfs , the following approximation is used. The function 

20/ - V bs ) = |[(Vi, + 20/ - V bs f' 2 - (20/ - y ti ) 3/2 ] (2.19) 


can ho approximated numerically in the range 0 < V ds < 20K and 0.6 < 20/ - V bi < 12.6 by 

/•’i ( ' 20/ - H,) = ^(20/ - V bs )V ds + — (2.20) 

7(20/ - Vi.) 

Where 


//(20/ - Vi.) = 1 - 


( 2 . 21 ) 


1.41 + 0.43(20/ - Vi,)' 

The detail of this approximation is discussed in sec. (2.6) and (2.7). The drain current in the 
linear region is given by 

*. r xr t r 

( 2 . 22 ) 

(2.23) 


j 1o f T jr ^ ^ds i T j- 

1(1 £ H gs ft 2 J 


Where 


Vr = 20/ + l/ /6 + 70 7(20/ - Vi.) 
and is defined at low V,/ A as shown, 



WITH,, W7r+ ? " 


Figure 2.2: - 



CHAPTER 2. THEORY OF MOSFET MODEL 


14 


1+ f(24>,-V b ,) [1 1.41 + 0.43(20, -Vi,) 1 ' (2 ' 24) 

In the model equations ”a” describes the lowering of drain current due to bulk charge in the 
substrate. The factor ”a” given by this equation (2.24) represents the bulk doping effect. 


2.4 Carrier Mobility Degradation 

Carrier mobility depends on effective normal field E^ and lateral field E L in the channel. As 
a first approximation, 


n Q 

M = f+aiijvi+^r (2 - 25) 

If Ei and E 2 are the electric fields normal to the surface at the Si— SiO^ interface and channel- 
space charge interface respectively, they may be determined from the following conditions: 


t s ito{E x E<T] — Qc and e s ,e 0 E 2 — Qg 


(2.26) 


where Qc and Q# are the channel and bulk charges, respectively. The effective normal 
electric field in the inversion laver will be 


\E n \ ~E x = 9* t Qc = ^E(v gs - 2 0/ - V fb - V) 






The lateral field El is given by 


\El\ = 


dV 

dY' 


from these equations, the mobility can be expressed as 




1 + 0[V gs - 20/ - V fb -V) + rtL 


dV 

dY 


Where 


and 


9 — a 


C a 




V =L 


(2.27) 


(2.28) 


(2.29) 


(2.30) 


(2.31) 


6 and 77 represent the electric field parameters for surface (normal) scattering and channel 
(lateral) scattering, of carriers in the inversion layer. In order to simplify this expression, it 
is assumed a constant channel electric field in the lateral direction, i.e., 

dY L ' 


(2.32) 
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This pcimits us to average the surface scattering contribution over the channel length as 
follows: 

* 1 1 /•Vrf 

I Jo ^ ~ Vfb ~ 2(j)f ~ V ^ dy - yr J 0 [Vgs -Vfb- 2 4 >f- V}dV 

= (V, gs ~ Vfb ~ 20/) — = V gs — Vt — + K\J~{2(pf — Vb s ) (2.33) 

With the preceding simplification, the mobility can be expressed as 

t* = — — — T (2.34) 

1 + Q[{V gs - V T f-) + K^f{24>j - V bs ] + r]V ds 

The general expression for the current in the linear region is then (y 0 = n a C ox ) modified as, 

-y W Vj 

h = -j-[V gs -V T - a-~]V d3 (2.35) 

where a is defined in equation (2.24) and. 

7 “ 1+0 ( -Vr- %+X y/Wj-V^Wi, 


2.5 Channel Length Modulation 


For short channel devices, the drain voltage corresponding to the onset of saturation for 

a given gate voltage is determined by the requirement that the lateral electric field at the 

pinch off point is at its critical value E c . For this field, carriers are at their limiting velocity. 

To find an expression for the saturation voltage V d3at , we equate the current calculated from 

the linear portion of the channel with the current at the pinch off point; hence, 

W Va , 

pC«-[V SJ -V T -a- 

= !tWC„[ V„ - V, b - 24/ - U,. t - Ii(V int + 24,, - V b ,C]E b 

= n WC„F 2 (V „ , 24, - V b ,)E c . (2.36) 

For usual values of voltage, channel length, and thin oxide thickness, we can write with 
acceptable accuracy. 


F 2 (V g „2ct> f -V b3 ) = 

V g3 — Vfb — 20/ — V dsat — K(V dsat + 20/ — Vbs ) 1 / 2 = V g3 — Vt — aV dsat 

W Va » 

nr llrv.. -Vr- ]V d ^ 


(2.3^ 
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from which we obtain the saturation voltage V dsat , 

v dsat = + E C L - + (E C L) 2 ) ( 2 . 38 ) 

a v a i 

For large values of channel length, E C L » V gs — Vr , and neglecting the bulk doping effect 
(a=l), we obtain, 

V dsat ~ F ss - F t . ( 2 . 39 ) 


The saturation current l dsat is determined by inserting Y dsat into equation (2.35) in place of 
F d4 


where, 


7 = 


j- 7^ fi/ f ^ Fisai i-j^ 

*dsat — r [*'gs Fr ® 2 j * dsat ■ 


7o 


1 + <?[(F, S - Ft - V) + A '\/(% - H s )] + ??F dsa , 

Channel length modulation for V ds > V dsat , results in the relationship 


( 2 . 40 ) 

( 2 . 41 ) 


I d = -~[Vg S - F r - a^}V dsat ( 2 . 42 ) 

Where derivation of L/ is given below A simplified treatment can be made using the follow- 
ing hypotheses. 

1) Mobile carriers have a saturation velocity v L at the pinchoff point in the channel. The 
electric field parallel to the surface at this point equals E c . 

2) Mobile carriers between the pinchoff point and the drain are pushed away from the surface 
by the normal component of the electric field and are collected by the drain diode at the 
depth Xj. These mobile carriers affect the space charge layer width in the same manner as 
in the same manner as in the BJT. 

3) For the sake of simplicity in the mathematical treatment, we assume that the carriers are ! 
spread uniformly in this region over a depth ranging from d for the pinchoff point to xj for 
the drain. 

Then the set of equations in the drain-channel space charge layer is 


d 2 v 

dy 2 


€ C 6 S 


■( qN D 


vl 




d 


( 2 . 43 ) ; 

( 2 . 44 ) 
( 2 - 45 ); 


X = 
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It is assumed, a voltage drop between source and pinchoff point equal to V dsa t, and boundary 
conditions 

v=() at V=i;, Sf „. and = E ( « % 
v=l at V=r f/ „ 

From the hypotheses and boundary conditions, the solution for the set of equations (2.43), 


L 




uyA 


j | 5 i 02 dj -* 

/ D / / /' Vr\<ie.£^-/ . 


K- JL- J 


j • ■ * M y-dy: a , 


V 0&> V DSS 


Si N 


Fit;. 1 . Section of the MOS transistor. 


(2.44) and (2.45) is 


r/iV ;) 


s 


mmmmmmwmmmMmmpmm gate 

- \5t 4 £n |h s '°2 

cw*i 


Spoc^-cKortj* 

BULK 


/*[! + 


iglire t 3 ■ ^ Normal and lateral fields. 
2 1 ds Xj 


.{log-A — 1 + — )] + lE e (2.46) 


2(„e ti ' L ~ ' qN D Zv L (xj - d) 2 '^ 3 d x 3 

This formula shows dearly the effect of drain-junction depth Xj on MOS transistor output 
impedance. If .r, increases, V,i s - Vdsat. decreases and there is decrease in output impedance. 
Output, impedance is inversely proportional to drain-junction depth and is maximum for 

xj = d. 

For the usual diffusion . :r 7 >> d and (2.46) gives 

2 /d., 


r \/ ( l N n 2 . 

Vtls ~~ VdHat — 7T } L 1 + 

<)£ si 


q N D Zv LXj {, °»J- l " + lE ' 


.v/:-' „ ._ and v c = e c L 


subsit ute . 1 H - ^.wn, 


0 = A/L?l^[l + Bids] 4" 2V c l— 2(Vd,, Yisat) 


•(2.47) 


(2.48) 


— 2V(. ± S(4Vc> + 8/1(1 + BIds){Fds Vdsal))-. 
1 := L( 2 / 1(1 + BlZ) 


(2.49) 



L L [VI + 2A(1 + Bhsadiy^ - Vtoat )] 1/2 - Vc 

^•(1 + Bldsat ) 


( 2 . 50 ) 


5 mobility modulation factor 9 factor consists of two parts. The first part, 6 0 , is caused by 
conventional mobility degradation due to surface scattering of carriers in the inversion 
;r. The second part, 9\jL. is interpreted to be the effect of intrinsic series resistance r t in 
source or the drain on the apparent mobility . The rj factor also consists of two parts. 
? rj 0 part is believed to exit due to the approximations used in deriving the mobility 


lendence on V ds , and its value is very close to zero. The second part, r/i/L, is due to 
nnel velocity saturation effects (channel scattering of carriers in the inversion layer). 



Ad • \J 


xiwo wa lutnusic source AJND/OR Drain Series 


Resistance On The Carrier Mobility 


rhe parameter Q\ is caused by the device ’’intrinsic” resistance r, that is due to current 
rowding in the drain-to-source and drain to channel regions, The results are valid for all 
•alues of V ds . Linear relationship (ForV^ << V g , - V T and V g , > V T ) 


h 


7 oW 

L 


\v,. - W,, 


here, 


id 


> simplify these equations, we get 


K = K - K> 



(2.51) 


(2.52) 

(2.53) 

(2.54) 


VL = Vds ~ IdTi. 


(2.55) 


in terms of V ds and V S3 , 

' 7o W 

Id = i + %(v gg - v T ) T ^ 5 ~ v ^ Vd3) ^ 2 - 56) 

% 

iere 6 X is defined as 


0\ = 70 Wri = constant. 70 . (2.57) 

:e ri is proportional to 1/W,this imply that 9 X is independent to W and 9 X is only 
endent of 70. o v < 



Intrinsic source and drain series resistances. 



Approximating Fi{V ds > 20/,- K^). 
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2.7 Numerical Approximation Of Function Fi(V ds , 2cj) f - 
' " V h8 ) 

The accurate approximation of function over a reasonable range of V ds and 2<p f - V^.Let 


V x = 2 d> f - V bt . 


For small V dx , 

a(u„ v x ) = |[(w, + v x )’i* - (iv) 3 / 2 ) 


by expanding it, 


Fi(V ds , V x ) = {V x ) l ' 2 V ds + 0.25 {V ds ) 2 /(y x fF + 


(2.58) 


(2.59) 


(2.60) 


This expansion is valid for only V ds « V x . The expansion is modified to 


Fi (V ds ,V x ) = (V x ) l / 2 V ds + 


0.25(V r ( f i ) 2 g(Vy) 

(V x )L2 


( 2 . 61 ) 


Where g(V x ) is determined by expanding this equation to give the best fit to Fi(V ds , V x ) over 
the desired voltage range. The parameter g for each fixed V x is determined the expression 
(V x )P 2 V ds + 0.25~^|i. can be accurately approximated by 


5 = 1 - 


1 

Pi + P 2 Vx 


(2.62) 


Where Pi and P 2 are determined by linear squares fitting over the range of V x . The result 
Is Pi = 1.408 andP 2 = 0.4333. 


2.8 Numerical Approximation of Function F 2 {V gs , 2<t>/ - 

V bs ) 


To find a good approximation of F 2 {V g3 , 2 <j>f-V bs ) over the usual' values of voltages .Assume, , 
F 2 = 2 (j) f - V k)l . Hence, 


1 F 2 (V gs , V x ) = V g3 - V fb - V d sat ~ K(V dsat + V x ) 1/2 - 2 Of 

= V gt - V fb - 2 4> f - V dsat - KyJVxll + ^] l/2 (2.63) 

Expanding the factor [1 + using Taylor series and taking first and second term only, 



5.8 Subthreshold Region 


Vhen gate voltage us below the threshold voltage and the semiconductor surface in weak 
iversion. ihe coi responding current is called the subthreshold current. In weak inversion, 
he drain current (dominated by diffusion), 


(in n(0) - n(L) 

Id = -qAD n — = qAD n ---..' ... \ L 

cly L 


(2.63) 


/here A is the cross section of the curernt flow, and n(0) and n(L) are the electron densities 
n the channel at. the source and the drain, respectively. These electron densities are given 

>y. 

n(0) = ripoeP** (2.64) 

n{L) = n vo e^^ v ° (2.65) 


chore is the surface potential at the source. The weak inversion surface field E s is given 
>A 

E,, = -Qb/cs = \j2qN A ip s /e s ■ (2.66) 

Subsituting Eqs. 2.64, 2.65 into 2.63 gives 

h = ^(f )^(^) 2 d - e-? v °y*-.mr 112 

Where q is the electronic charge; //.„ is the effective electron mobility; L is the channel length; 
vV is the transistor width; The Schematic diagram of a MOSFET is shown below, 

Because the current in a MOSFET is transported predominantly by carriers of one polarity 
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only (e.g., electrons in an n-channel device), the MOSFET is usually referred to as a unipo- 
lar device. The MOSFET is a member of the family of field effect transistors. Although 
MOSFETs have been made with various semiconductors , the most important system is the 
Si — Si0 2 combination. We first consider the basic device charactristics of the so called long 
channel MOSFET; that is, the channel length L is much longer than the sum of the source 
and drain depletion layer widths ( W s + Wd)- This serves as a foundation to understand 
short channel, that is,L <= (W$ + W D ) , and related MOSFET devices. For short channel 
devices, the full effect of Q B on the threshold voltage is reduced, because near the source and 
drain ends of the channel, some field lines originating from the source or drain terminate at 
the bulk charges in the channel region for V D = 0. For V D > 0,the depletion region near the 
drain expands further . Note that the horizontal depletion layer widths y s and y D are smaller 
than the vertical depletion layer widths Ws and W B , respec tively, because the transeverse 
field strongly influences the potential distri bution at the surface. 

Because of the reduction of the bulk charge Qb, surface potential for a given gate voltage 
increases, leading to an increase of subthreshold current . The surface potential can be found 
from the following expression: 

V G ~ V FB = i’ s + ±-f(qe.N A {t/>. + N A (ij>, + V BS )/2))(l + (2.68) 

Ci L — y D — y s 

where W D and W s are given, 

H d = 7(-%-(V d + + F B s)) pm (2.69) 

v qNa 

IF.s- can be get by putting Vp=0, 

Ws = ^Wa (Vk + Fbs)) (2 ' 70) 

ys = i(~(V D + V k + V BS ) (2.71) 

vd = - *)) < 2 - 72 ) 

The subthreshold current is given by 

I D = iJ. 0 {- r — (^ L ) 2 ( 1 “ ^P (-PV D ) exp(^)(^) 1/2 (2.73) 

L - y s - Vd Ap) ■ /v a 

Equation (2.73) is identical to equation (2.67), except channel length is replaced by the ef- 
fective channel length: 


L e jf = L - y s - yD- 


(2.74) 



Chapter 3 


The MOSFET Model 


3.1 The Model 

The model equations are derived - — - in chapter 2. 

3.1.1 Saturation Voltage-V^ so< 

The saturation voltage given by eq. 3.1 below is used to determine whether a device is in 
the linear or saturated region. The equation is derived in chapter 2 (eq.(2.38)) 

' = v " ~ Vr + ) 2 + '-a (3.i) 

a v a 

Where 

\ ’ gs = gate-to-source voltage. 

\ r ds— drain-to-source voltage. 

Vbs= bulk (or substrate)-to-source voltage. 

V T = V t (L, W. Ti ds) is the device threshold voltage. For small value of V ds can be (<100mv), 
V T can be determined. 


V e = E C L = —L (3.2) 

A^o 

= — (3-3) 

t-O^OX 

E c critical field in the channel (for this field, carriers are in their limiting velocity v L ) in 
V/cm. 

carrier mobility at low gate and drain fields in cm 2 /V.s. 


23 
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channel length (wafer) in cm. 

/ channel width (wafer) in cm. 

0 gain factor in U/V. 
tx oxide thickness in cm. 

, permittivity ol free space in F/cm. 
„ relative permittivity of SiOo. 


,, relative permittivity of Si; and 

'his ('(jiiation 3.-1 has been derived in chapter 2(eq.2.24). 


(I — 1 "h 


0.5 A' 


.fl t 1 

- vy 1-41 + 0.43(20/ - v b ,y 


(3.4) 


diere 


1 slope 1 of Vr versus — V bs ) curve. Equation for K is derived in chapt< . 2. As is 

vident from eq>(2.23) 

T — 2 <Xf + \ fi, + A — H«)- 

dins, 


K 


V T2 - V T1 


- V bs2 ) - yf(2<p/ - Hsl) 


(3/5) 



-ri, Vn , \ /mi, and are as shown in lig.(3.I). 

Upj twin' the Fermi potential level in the ’’bulk” 2 (j>f — 0.6V for bulk doping ol 1.2cl5 cni 3. 
rhe value of K is dependent on the doping profile. ’For uniform doping, in more practical 
use wlu'ii threshold tailoring results in nonuniform doping, K may be treated as a constant 
,o a good approximation over the voltage 'range of interest (Vd s < — l.o and V r i s > 0V ). 
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3.1.2 Linear Region 

0< V ds < V dsat and V gs >= V T . From chapter 2(eq.2.35) 

T t 7 Ids it r 

I d — —j^Wgs -V T - a — ]Ku 

where, equation for 7 has been taken from chapter 2 (eq.(2.41)). 


7 


Equation for 77 and (3 , 


lo 


1 + 6{V 9 , - y r - \ + Kj(2<t> f - V bs ) + rjV ds 


e = e 0 + 


Oi 

V 


, m 

ri = r ) 0 + —- 


3.1.3 Saturation Region 

i <is 7 1 dsat and \‘gs 7 = ty 


7^ fj r 7 / 

- lr - a-j-]Vdsat 


where 


(3.6) 


(3.7) 


(3.8) 

(3.9) 


(3.10) 


r r r [K 2 + 2A(1 + BI dsat )(V ds - F dsat )] 1/2 - V c 

Jj r r=r JU Ju , 

>1(1 + Bldsat) 

(3.11) 

a^ nl2 (V). 

^o^si 

(3.12) 

B = 2( '^- 1) (A)-l. 

qNv L WDj v 

(3.13) 

q=electron charge. 

N=impurity density at the drain space charge region. 

Injunction depth(cm). 

d=mean inversion layer thickness (cm). 


Vl r 
v c = — L. 

Vo 

(3.14) 

lo 

! ~ 1 4- e[(v g , -k T -V) + *V( 2 */ - 

. (3.15) 


7 


(3.15) 
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(3.16) 

(3.17) 

3.1.4 Cut off Region 

Vp S < V T . (3.18) 

Id = 0. (3.19) 

The model parameters 7 0 , K, 6 and 77 obey certain relationships with physical parameters as 
are explained in the following. 7 0 is the gain factor and is related to the physical parameters 
by the relationship 

7 O = 1-^oOox (3.20) 

where 

C ox is the oxide capacitance per unit area in F/cm 2 . . 

It is emphasized that the mobility modulation factors 6 and 77 modify 7 according to model 
equation (3.15). The factor ”a” describes the lowering of the drain current due to bulk 
charge in the substrate according to eq.(3.10). The effect of current modulation due to the 
short and narrow channel effects is embedded in the model parameter K. Note that when a 
is set to 1 and 0 and 77 are set to 0 , 

Id = ^{V gs - V T - '-f)V d$ forO <= V ds < V dsat (3.21) 

and 

Id = - V r ] 2 forV ds >= V dsaU (3.22) 

Which are the classical equations for a simple MOSFET model. 

3.2 For Depletion Mode 

The turn on characteristics, l ds versus V gs measured with a small drain voltage, are shown 
here for enhancement and depletion devices with the same physical dimension. A threshold 
value V T of gate voltage is defined by matching the experimental curve to a straight line 
segment asymptotically at the region of maximum slope, extrapolated to zero current, in- 
tersecting the axis at V T + V ds /2, from which V T is determined. When the gate voltage is 


6 — 0 o + 


V = Vo + 


Li' 

Hi 

Li' 
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below V'r Kn, the device is in pinchoff, and the current varies as [V g> - V';-) 2 according to 
the simple theory. The depletion device behaves somewhat differently in this "soft turn-on 
region" and has a definite (but small) current level at V g , = Vp due to substsface conduc- 
tion in the implanted n-t.vpe layer. The straight line segment matches the depletion device 
current very well at V gs — 0, which is the usual bias condition for e. depletion type load 
device in standard logic circuits, so this threshold definition is meaningful for modelling such 
circuits. The slopes of the straight line segment are similar for enhancement and depletion 
mode devices. 



Turn-on .'harsctcrUtici ol f;MD and DM D device* nhowinji 
threshold definitions. 


Figure 3. II 

The current behaves of the depletion mode devices with shallow implanted channels can 
be adequately described by 3.1-3.17 equations, which represent an enhancement surface 
channel device when defining the depletion device threshold Vp as explained above. The 
turn-on characteristics, Ij,vsV g , measured with a small drain voltage, are shown in fig.(3.il). 


3.3 Extracting Model Parameters 

A parameter extraction method for the above model is described below. It uses conventional 
methods to extract Vp and K, and uses a linear technique to extract 7o, 0 o .0 u r] O) and r? i- 
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The measured data were taken from reference (4) in Bibliography with bias voltage of VJ, S 
from -1.5 to 8 V, V gs from 0.5 to 4 V above \V T \, and V ds from 0 to 5 V. For enhancement 
'mode devices with channel lengths from 1.0 to 13 fim, channel widths from 2 to 27 flm, 
substrate doping = 1.2el5cm 3 oxide thickness = 45nm, and junction depth = 0.5 fim. The 
saturation velocity vj_, for electrons was assumed = 15e6cm/s and the mean inversion layer 
thickness d was assumed to be 10 nm. For small values of V ds (= 0.1V) and large values of 
VgstVgs - V T ) » I</{2ct> f ~ V bs ), 


L 1 + 6{V gs - V T - fy) ds 


This also can be reduces in this form, 


(3.23) 


Rds 


V* 

Id 




Po"V gs -V T 




+ 6 ) 


(3.24) 


where 



(3.25) 


(The programme for MOSFET is presented in sec.C.l of Appendix C.) 
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3.4 Results 

The MOSFET model for short channel is described above. The drain current Ii>5 is ploted 
against V* for different regions or different conditions which are dependent on the values of 
V 9 s and V r . Fig. 3.1. 3.2. 3.3. 3.4 shows the drain current characteristics using this model. 

Fig-3.1 (a) ‘and (b) shows comparison of drain current and drain voltage characteristics in 
enhancement mode device for constant 0 O and different channel widths. 

Fig 3.2(a) and (b) shows comparison of drain current and voltage characteristics in enhance- 
ment mode device for constant channel length, channel width and different 6 0 . 

Fig. 3. 3(a) and (b) shows comparison of drain current and voltage characteristics in depletion 
mode derice for constant 6 0 . channel width and different channel lengths. 

Fig. 3. 4(a) and (b) shows comparison of drain current and voltage characteristics in depletion 
mode derice for constant 6 0 and different channel lengths and widths. 

The slope of the line ^ versus [ — ~ — ^ j represents ^ and y-intercept shows This is 
shown in fig. 3.5 (a). Family of these curves give values ^ and A. 

To plot ■— versus L m , it is found that slope of this line is equal to and the x-intercept 
is equal to A L. where A I = Design channel length -Actual channel length . This is shown 
in fig. 3. 5(b) for neighbouring transistors with constant W and different L. 

AH’ is determined by the measure of x-intercept of the line drawn between and W m . This 
is shown in fig. 3. 6(a). Fig.3.6(a) is ( 3 q vs W with x-intercept, where AH'= design channel 
wddth.. -actual channel width. . 

Fig 3.6(b) is a plot of the line 7^ versus [— -A-p^- ] for constant length and different widths. 

The dope of these lines is ± mi y . intw " 1 

A> 

The e„ is determined as the slope of the line which is drawn between £ and J- and y- 
intercept is equal to r, + R,. This is shown in fig.3.7(a). * ^ 

Fig.3.7(b) is a plot of r, versus & and slope of this plot gives the value of 
Information obtained from fig.3.5(a) is used in fig.3.7(b) i.e. A and A where 1 vs - is 
plotted in fig.3.7(a) the y-intercept gives r. + R, where R. is the parasitic resist^ce nTthe 
source or the drain and 6 0 is the slope of the line. Plot of v vs ft shown in fig. 3.7(b) is a 
straight fine with extrapolated intercept gives t*. tj, is obtained from the slope. 

There is good agreement between measured and calculated data. In this model, a parameter 
extraction method has been used to extract 7„, 0 O >% and r/j. It has been found that the ac 
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curate of the simple MOSFET current model car be significantly improved by including the 
effects of the entire bulk doping term, the dependence of carrier mobility on gate voltage, 
drain voltage, intrinsic source and drain series resistance and channel length modulation. 
Also shown is detailed derivation of the model equations and a method of extracting model 
parameter. 
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EMD 

W/L=2.16/13.4 



Ids vs Yds 


=.008 

1^=25 



Yds volts 


(a) 


D j =0.5 ji m 
n = 1.2 X itf^cm 


I 


V tl = 80 m* 

dsat =lmV 

V = 40m V 

t2 

V = 45 mV 
T 


\s =' 3Volts 



v j = 15 X 10 6 cm /s 

d = 10 nm 

t = 45nm 
ox 

<b = 0.3 V 

f 

^ s2 = -8.0V 


M,sl = - 1 ' 5 V 


(b) 


Figure 3.1: (a.)-{b)I d vsV ds characteristics of MOSFET for different widths 
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Figure 3.2: (a)-(b )IdVsV ds char acteristics of MOSFET for different parameter B 0 
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Figure 3.3: (a)-(b )I d vsV ds char acteristics of MOSFET for different lengths 
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Figure 3.4: (a)-(b )I d vsV ds char acteristics of MOSFET for different widths and lengths 



























Chapter 4 

'he Model Of The Sub-Threshold 
Characteristics Of Poly-Si TFT 

1 The Model 

; electronic properties of a-Si and their relationship with deposition conditions are rela- 
ly well characterized owing to the huge amount, of fundamental and applied work doin' 
this material in past years. Models correlating TFT performance to the gap density of 
es (DOS) have been formulated. 'Most of the characteristics work on this material has 
i done for applications in integrated circuits such as a resistors and interconnections for 
eh a detailed knowledge of the electronic properties of the material is not crucial. The 
'nnination of the density of states of poly-Si, based on the analysis of the transistor 
■aetenstics. Density of states is continuous in energy and depends strongly upon the 
morphology. Based on these, a new model for the characteristics of poly-Si TFTs is 
.iced. We assume here negligible free carrier concentration and small drain voltage V D 



13 

pit; ,J { Density of states lor three samples' ( I ) poly nlaon thickness 0 hjt m 

( A ); (2) polysilicon thickness i 5//m ( D). (3) sample A niter b)diogena 
tion (C). 
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and the polysilicon film consists of 2-dimensional grain of length a. Poisson’s equation near 
a grain boundary at y=0 as shown in fig.B.l of Appendix B; 

cfrp _ S(x)8(y) 

dx 2 dy 2 e s 1 j 

where x is the direction perpendicular to the channel, ip is the electrostatic potential , 
S(x) is the grain boundary charge per unit area, 8(y) is Dirac’s 8— function and e s is the 
semiconductor dielectric constant. As for the boundary conditions for V D = 0, symmetry 
considerations require: 

dip . 

~dy^ ±a/2 = ° 

Thus, for small Vu it is reasonable to write: 

-Ys. 

dy \ ±a/2 L ^ 4 ' 3 ^ 

TFT subthreshold characteristics i.e. V D ~ 0 or V d is small. Where L is the channel length. 
By integrating Poisson’s equation with these two conditions from to |: 


where 


_ 1 W 2 

4’ = - rp{x, y)dy 

a J —a/2 


and e(x)=S(x)/a. When the Debye length Ld is much longer than the grain size, ip = 
ip(x, 0) = ip{x . + a/2) and the problem is reduced to a 1-dimensional one, where e(x) is 
the average charge density per unit volume. Since the polysilicon used for TFTs is either 
undoped or lightly doped, and a ~ lOO-lOOOA. From eq(4.4) 

0 &P &ip dip e(x) 

dx ; dx 2 dx e, ' 


d_M 2 

dx K dx 


dil’ e(x) 
dx e. 


Integrating above eq. 


r = r 2^dA ix 

Jo dx dx Jo dx e s 

ii — -) 2 =f/>w 

^‘° f =zL e( ' p)i ' p 


(4.10) 
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In n-channel TFT, for the DOS in the upper*of the gap, it can be assumed without loss of 
generality, and consistent with the experimental data, taken from reference (7) and (8) of 
Bibliography: 

N(E) = N G exp((E - E C )/KT G ) (4.11) 


where N G is the density of the gap states immediately below the conduction band edge E c 
and I\T g is the characteristics energy of the exponential band tail. By multiplying both 
sides of eq.(4.4) by dip and integrating in x from the interface (x=0, il> = V ; «) to the back 
surface (x=d, -0 = 0) We get: 




(4.12) 


To approximate 
temperatures. 


the Fermi function to a step function, which is reasonable for low enough 




= F 2 L(ii,) 


(4.13) 


where: 


L(i> s ) = exp (qip s /KT G ) - qip s /KT G - 1 (4.14) 

F 2 = 2 N T KT G /e s (4.15) 

N t = N g KT g exp ({E F - Ec)/KT g ) (4.16) 

The channel sheet conductance (G) of the device is, 

G = G 0 + Go/d [ d (exp (qiP(x)/KT) - 1 )dx = G 0 + G 0 /d (*' ~ 1 dij , (4.17) 

Jo Jo dyj/dx 

where Go is the flat band sheet conductance. We assume here for simplicity n-type conduc- 
tivity in the flat band conditions: 


Go = qpdN c exp((E f - E e )/KT) 


(4.18) 


Eq(4.17) by differentiation gives: 


dG _ -Go exp {gxjjs/KT) - 1 
dip s d (dip/dx)x — 0 

knowing, 

V G ~ Vfb — h'ox + fps 


FVom eq.(^) 2 |* =0 = F 2 L(ip s ) Assume, 


(4.19) 


(4.20) 



CHAPTER 4. THE MODEL OF THE SUB-THRESHOLD CHARACTEREISTICS OF POLY-SI TFT41 


(^) 2 |*=° = F 2 exp(qip s /kT G ) 

(4.22) 


dG Go exp(qip,/kT) 

dip, d Fcxp(^) 

(4.23) 


dG G 0 , . qip, . 

a p^t.IkT 2kTa ) 

(4.24) 


2|-D) 

(4.25) 

neglecting the effect of interface states, 



- V FB = tox 7~ x=0 + 

e 0 x a# 

(4.26) 

For small -ip s 

V V ~i £ ° d ^\ 

VG — V FB — “ |x-0 

&OX 

(4.27) 


V G V F b ~ F exp 

€ox j&KA q 

(4.28) 


qips v G — V'/.’/i Vg — Vpb 

eiP 2 kT a ~ t 0x ‘f o F ~ Vo 

(4.29) 

using above equations we 

have 



— = gsf( V g ~ Vfb )^~ 1] 

dip, V 0 

(4.30) 

further 

dV G __ t ox e s Fq # a 

dip, e 0 2kT a FK 2kTc: 

(4.3i) ; 


dV G t 0X e s Fq V G - V> B 

dip, e 0 2kT a 1 V 0 

(4.32) 


d0 s _ e 0 2kT G , Vo , 

^V; t ox e s Fq V G — V FB 

(4.33) j 


dG dG dip, 
dV G ~ dip, dV G 

(4.34) 


dG Gotox^kTo ,V G — Vp U yWa_ 2 ) 

dV G qdF' 2 e s t 0X V 0 

1 

(4.35) 

1 

using 

F l = 2 NrkTcps 

(4.3G) 


dG GqCox / bf; ~ Vpu y'£Ta_2) 

dV G qN T dt ox V 0 

(4.37) 

[ 
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Integrating 

r G 0 e OI 1 (Vg ~ Vfb ) 2Tc / T ~ 1 

° ~ qN T dt ox 2 T G /T - 1 (Vo) 2r °/ r-2 

(4.38) 

r c Got,, T (V 0 - V FB ) 2T °' T -' 

° qN T dt or 27c - T (V 0 y r o/r -2 

(4.39) 

F 2 = 2N t KTg 

(4.40) 

( Go = qndN c exp((Ep - E c )/kT) 

(4.41) 

Nt — A’c kT G exp(( Ep — E G )/kT G ) 

(4.42) 

Y _ £s*ox ( 2 NrkTc y / 2 

, ^OX 

(4.43) 

y 2^-2 _ Y^-- 2 ^ T kT G 

€ox €3 

’(4.44) 

Gq qfidA'c r Ep — Ec Ep — Ec, 
Nr = N c kT c eZpl kT kT a 1 

(4.45) 

G 0 _ qpdNc ,E f - E c ,T C 
, N t N G kT G P[ kT G [ T )l 

(4.46) 

G 0 qpdNc , Np 

• N T ~ N G kT G [ N G kT G ) 

(4.47) 

K _ G 0 e ox , T 1 2ia_ 2 

0 qN T dt OI { 2 T G -T ){ Vo } ' 

(4.48) 

K e 0 X N r qpdN C( N t T f e OI e, 1 ,^_i 

qdt 0 X N G kT G K N G kT G } 2 T c -T K e s tJ K N T ’ K 2kT G ’ _ 

. (4.49) 

’ R 0 -qp Nc ( 2 T • G _T^N G q 2> ( 2e i t 0 X 2 k 2 T G 2 N G ) * ? 

(4.50) 

where 


Vo = C)t 0 X F/€ 0 x 

(4.51) 
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It is above derived that the integration of eq.(4 3T) gives: 

G — G? — Rq(Vq — Vi?b)^ 2Tg ^ t ~^ (4.52) 

where A' 0 - ^7V t (T/(2T G - r))(2e,/^?^ (The pro- 

gramme for Sub Threshold characteristics of TFT is presented in sec C 2 of Appendix C.) 



CHAPTER 4. THE MODEL OF THE SUB-THRESHOLD CHARACTEREISTICS OF POLY-SI TFT44 


4.2 Results 

The subthreshold characteristics of polysilicon TFTs are explained by neglecting the discrete 
grain structure and assuming a spatially distributed concentration of defects. The mathe- 
matical modelling of the dG/dVa vs Vq — Vfb equations presented here is more direct and 
gives a better idea of the approximations involved. At different temperatures, plots gives 
the validity of the approach. 

A plot of log(^) vs log(Vc - Vfb) should give a straight line of slope 2^ — 2 and the 
characteristics parameter of the exponential band tail To can be determined. This is shown 
in fig.4.1. 

A plot of LOG(G) vs Vg is shown in fig.4.2 which gives G-Vq characteristics for T=310 K. 



LOG(G) 
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Figure 4.1: 
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Figure 4.2: 





Chapter 5 


Modelling for accumulation mode of 
Poly-Si TFT 


5.1 The Model 

A schematic of polysilicon TFT is shown in fig. (5. 1). 

Polysilicon thin film transistor consists of polvsilicon film usually grown by low-power chem- 
ical vapour deposition (LPCVD) approximately lOOOAthick Gate oxide (1000A) and polysil- 
icon gate electrode (3000A) is deposited on the polysilicon film. The source and drain are 
formed by doping using ion-implatation and subsequent annealing. Polysilicon TFT is driven 


gore $,0 2 gate electrode 

S -V d 




It : 


source 

electrode 




V(x) 

• 5 — 

- L 

| i 

drar 

- n + j electrode 

i 1 

L*J 


polycrv" t alline silicon inversion 

film layer 


Cross section of polycrystallme thin-PJm transistor showing the 
symbols used in the text 

fio] 

fig. (5.1) 

into accumulation mode when gate voltage V g exceeding the threshold voltage V th is applied 
to the poly-Si TFT. The carriers, which are electrons in the n-channel mode, are induced 
near the polysilicon film surface, and a strong inversion layer is formed. The carrier concen- 
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tration n(x) in the strong-inversion layer for potential V(x) at the distance x from the source 


electrode is 

n(x) = C ox (V g - V th - V(x))/q.d (5.1) 


where C ox is gate oxide capacitance, and q and d are the electronic charge and inversion layer 
thickness, respectively. For strong inversion conditions within the grains, electron trapping at 
localized grain-boundary states produces potential barriers that effect the electron transport 
along the channel. The barrier formation is similar to that at grain boundaries in bulk 
polysilicon, except in the channel ip(x) is influnced by V g as described by the two dimensional 
form of poisson’s equation. 

Consider the potential variation near a grain boundary in the channel as shown below. 

It is assumed that away from the grain boundary (y > yj) the electric field is vertical (in 



Cross section of effective inversion layer showing typical grain 
and grain boundaries, which are assumed to be perpendicular to the 
electron flow. 

T ,. r -19 
Figure 5 II . 

the x-direction), and the true electron distribution n(x) is well approximated by n(x) over 
the effective inversion layer, 0 < x < x i(e// ). In this region poisson’s equation simplifies to 


&)\, = —(«(*) + N a ) (5.2) 

ox* e s 

where ip is electrostatic potential. It is assumed that that trapping nearly depletes the region 
(II) of free electrons;hence 


( 


d 2 ip, 


dx 2 ' 111 ' K dy 2 • e s 
Assuming that, analogus to the gradual channel approximation, the trapped electrons at the 
grain boundary typically create only a small perturbation on the x-coijiponent of electric 
field, so 


+ (S)I 


II 


-(AW 


(5.3) 
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where partial derivatives are evaluated anywhere in the regions indicated. Prom the comi- 
nation of (5.2)-(5.4), the corresponding approximation for ip{x) 

(fjoi// ~ (“M*) ( 5 - 5 ) 

with the boundary conditions 

ip(x,y = y d ) = ipj(x) (5.6) 


and 


dip 

dy 


y=yd 


= 0 


(S.7) 


The solution to (5.5)-(5.7) is 

ip(x, y)\n ~ (y - y d f + ipi(x) (5.8) 

and hence grain boundary potential 

ip(x) = ipi{ x) - ip(x, 0)1// ~ T}^)y± (5.9) 

ze s 

By the conservation of charge in the vicinity of the grain boundary 


Qgb = -2 qn{x)y d (5.10) 

The trapped charge density depends on the distribution in the energy gap of localized grain 
boundary states. By approximating this distribution by a delta function yielding N st states 
(traps) per unit area at energy levels Ep. Then eq (5.11) 


Qgb = 


~QN s t 

1 + larp ^-gr ' - } 


(5.11) 


where Ep is the fermi level and factor 1/2 reflects the degeneracy of the localized states. 
The position of Ep relative to Ep is defined by ip(x) and the electron density in region I, 
n(x) 


[Et ~ Ep\y = 0 


[E t — Ei] + qip(x) — kTln( 


n(x)_ 

Tli 


(5.12) 


where E , is the intrinsic Fermi level and n; is the intrinsic carrier density in silicon. V g is 
sufficiently high, ip{x) decreases with increasing V g , and hence Qgb — ~qN st is independent 
of V g . From eq (5.9) and (5.10) 


ip{x) ~ 


q_Nj 

Re un ( t) 


(5.13) 
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le induced carrier are trapped by trap levels at the grain boundary, which form a potential 
rrier. When the gate voltage is high, all carrier trap levels are filled. The expression for 
rrier potential height 

ip(x) = qN st 2 /8qi(x) (5-14) 

h ... ' In the other case, at low gate voltage, the carrier 
ip levels at the grain boundary are not completely filled due to the low carrier concentration 
iuced by the gate voltage. Therefore Fermi level is near the position of the carrier trap 
rel. The potential barrier height 4>(x) is approximated over a wide range of gate voltage 


ip(x) 1 = 1 + (qN st 1 /8^n(x)) 1 (5.15) 

us results in a virtual fixing of the potential barrier at V’o which is ~ 0.55V corresponding 
half the energy gap of silicon. Assuming that the carrier concentration taking part in 
3 conductivity is not the total grain-carrier concentration n(x), but an activated carrier 
ncentration n(x)’ defined by. 

n(x)' = n(x) exp (—qip(x)/kT). (5.16) 


r the interpretation of the hall effect it was assumed that the carrier concentration taking 
rt in the conductivity is the total grain carrier concentration n(x) and the mobility is 
tivated and has the form p = pQe~ q ^ x ^ kT . 

hen the mean free path in the grain is short compared to the grain size, account should be 
-cen of the scattering in the grain. In the general case the total resistivity p has contributions 
>m both the grains pa and their boundaries p# . Then it is written, 


P = Pg + Pb ~ 


+ 


(5.17) 


qpG^o qpbn 0 exp(-q'ip(x)/kT) 
lere pc is the grain mobility. Now, it can be written as p{x)~ l = p g ~ 1 +{pb exp(— qij)[x) / kT))~ 
irrier scattering takes place at the grain boundary to a greater extent. Then in the grain, 
ability pb at the grain boundary is less than the grain mobility p g in good quality grains. 


p(x) 1 = p g 1 + ( p b exp(-qip(x)/kT )) 


-1 


~ ( p b exp(-qip(x)/kT )) 1 (5.18) 


le gradual channel approach gives the current- voltage relations for poly-Si TFT’s with the 
annel length L and width W in the region of (V g — VA) > Vj 

I = ( W/L)C ox [ d Pb{V g - V th -V(x)).exp(-q^(x)/kT)dV{x) 

Jo 


(5.19) 
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where V d is the drain voltage. For the region of V g — V th » V d , the current voltage charac- 
teristics is defined as 

/ = (W/L)C oxi i*((V g - Vth)V d - V//2) (5.20) 

where p* is effective mobility and this is given 

p* — /j,b exp(— qip'/kT) (5-21) 

V'' indicates the effective potential barrier height in the region and is given by 

U' = a/(V g - V th + a/rk) (5-22) 

a = q 2 N st 2 d/8eC ox (5.23) 


(The programme for accumulation mode of TFT is presented in sec.C.l of Appendix C.) 
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5.2 Results 

r ig.5.1 shows changes of effective mobility as a function of gate voltage using equation (5. 21:. 
The effective mobility increases with an increase in the gate voltage. 

The I-Y characteristics were calculated using equation 5.21 for different gate voltage. Fig. 5. 2 
jhows I-Y charct eristics of polvcrystalline silicon TFT . 

Fig.5.3 shows changes of potential barrier height as a function of gate voltage using equation 
5.22. These all figures showa good agreement with experiment data (taken from reference 

;io)). 
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Figure 5.3: 





Chapter 6 


Conclusion 


The present work shows that the accuracy of the simple MOSFET current model can be sig- 
nificantly improved by including the effects of the entire bulk doping term, the dependence 
of carrier mobility on gate voltage, drain voltage, intrinsic resistance and channel length 
modulation. 

The MOSFET model presented here describes the drain current-voltage characteristics of 
surface channel enhancement mode devices in strong inversion and ion-implanted depletion 
mode buried channel devices. This model also presents a method of extracting model pa- 
rameters. In this MOSFET model, there is not included subthreshold conduction (soft turn 
on region) for depletion mode devices. 

In this thesis, the TFT model presented here explains the subthreshold characteristics and 
accumulation mode of polycrystalline silicon thin film transistors by neglecting the discrete 
grain structure and assuming a spatially distributed concentration of defects. To account 
for the observed I G — V G behaviour, it is essential that the density of states be continuous 
in energy. There is good fit of data with a DOS decreasing exponentially from the band 
edge towards the gap centre and with characteristics parameters, where expressions for G 
versus V G is deduced . The mathematical derivation of the ^ versus V G — V FB equations 
is presented, here. That gives a better idea of the approximations involved. The mobility at 
the grain boundary and the carrier trap density could be easily obtained, the mobility and 
trap density exert a strong influence on polysilicon thin film transistors characteristics. In 
addition, it is found from calculations that the effective mobility increases with increase in 
gate voltage. 

The following conclusions can be derived from present study: 
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the characteristics of the T.F.T. is similiar to that of the conventional FET, in so far as 
changes in the channel conductivity as a function of gate potential are concerned. 



Appendix A 


[’lie analysis ol MOSFET theory is divided into two sections: the first section gives a qualita- 
ive desci lption of the workings of the MOS; the second section concentrates on a quantitative 
liscussion which includes derivation derives the charactristic equations describing the MOS. 
Pho appioximat.ions made are used described below: 

1. Mobility of current carriers in the channel is constant. 

!. The variation of the channel thickness is small along the length of the channel. 

>. The thickness of the dielectric over the channel region is assumed to be much greater 
ban the channel thickness. 

I. Parasitic resistances (such as in the source) are assumed to be so small as to be negligible. 
). The channel is completely shielded from the drain, so no drain-to-channel feedback exists, 
i Doping of the 1 substrate is uniform and nondegenerate. 

The drain current consists only of channel current. Leakage currents are neglected, 
h The gat.e dielectric is considered to be a perfect insulator. 

). Extrinsic conditions which affect the conduction properties-such as oxide traps, silicon 
mrface states, interface energy states, ionic centers within the oxide, and work function dif- 
eieiices will be lumped together into a single effective charge term, Q 3S . Furthermore, Q ss 
s assumed to be constant, and located at the silicon-oxide interface. 


A.l Short Channel Effects 

> 

As the channel length is reduced, departure from long-channel behaviour, may occur. These 
departures, the short, channel effects, arise as a result of a two dimensinol potential distribu- 
tion and high electric fields in the channel region. When the channel length is reduced, the 


55 



APPENDIX A. 


56 


tion and high electric fields in the channel region. When the channel length is reduced, the 
depletion layer widths of the source and drain junctions become comparable to the channel 
length. On account of this, potential distribution becomes two dimensional because it de- 
pends upon two fields transverse and longitudinal fields. This results in degradation of the 
subthreshold behaviour, dependence of the threshold voltage on the channel length and bi- 
asing voltages, and failure of current saturation due to punch through. The channel mobility 
is a function of field due to increase of electric field. When electric field is much higher than 
carrier multiplication near the drain occurs. High field also cause hot-carrier injection into 
the oxide. 


A. 2 Qualitative Analysis 

Even though the following analysis refers specifically to a P-channel MOS on an N-type 
substrate (for convenience), the resulting equations are applicable to both N-channel and 
P-channel devices. There are three distinct conditions or regions occuring within the semi- 
conductor at the surface, that are important to MOS operation. They are accumulation, 
depletion, and inversion regions and are controlled by the external bias on the gate electrode. 
Generally, for an oxide passivated surface, surface states or energy states at the silicon-oxide 
interface act as ionized doners whose effect is the same as a positive applied gate voltage. 
Figure A.l shows the MOS structure used in this discussion. Drain-to-source voltage is as- 
sumed to be so small as to be negligible. The energy band diagram of a P-channel device 
under zero applied voltage is shown in figure. Here and in the following band-diagram , the 
intrinsic level (P = N = n;) is designated as E, and is assumed to be halfway between the 
conduction band energy E c and the valence-band energy E v . Because of the positive surface- 
state charge, negative electrons from within the N-type bulk are attracted to and accumulate 
at the surface (x=0). Accumulation results in a downward bending of the conduction and 
valence bands. 

The closer E c is bent toward the Fermi level, which is set by the substrate doping , the 
heavier the surface concentration of electrons becomes. Figure illustrates the charge density 
distribution. The positive charge per unit surface area ( Q ss ) must be balancd by the negative 
charge accumulated near The silicon surface (Qa)- (Charge distributions are approximated 
by the delta functions) If a small positive bias is now applied to the gate, additional band 
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Figure A. L: 
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bending and accumulation result. Again, the total positive charge must equal the total neg- 
ative charge so as to maintain charge neutrality (Q g + Qss + Qa — 0) - If a negative voltage 
is applied to the gate such that it just counters the effect of Qss, then no bending of the 
bands exists, a condition which is known as the fiat band case (bs = 0). Further application 
of negative gate voltage repels from the channel region the mobile electrons associated with 
donar centers, causing a depletion region to form. The band diagram and charge picture for 
this case are shown in figure. When electron is removed from its donar atom is left with a 
net positive charge. Thus the charge in the depletion region is shown as a positive charge, 
Q d . in figure. When the intrinsic Fermi level E, is bent just enough to intersect the Fermi 
level at x=0, the surface has gone from its initial N-type concentration to intrinsic (where 
P=N). (<p s is now equal to <pf.) Additional negative gate bias does notextend the depletion 
region so much as itinduces positive mobie holes at the surface. It is important to note that 
depletion region charge and mobile surface (channel) charge of the same polarity (positive) 
and thus their effects add. The sum of these two charge regions must balance the net charge 
stored in the oxide (Qss) and on the gate (Qg) so as to maintain an electrically neutral sys- 
tem. As the gate bias is increased in the negative direction, a larger and larger percentage 
of the charge within the semiconductor is contributed by the mobile holes. Until Ei = Ef, 
mobile electrons still outnumber mobile holes. Beyond this point , electrons are suppressed 
below the intrinsic level while holes in the channel region raised above this level. As will be 
seen later, conduction between the P-type source and drain (by P-type carriers) does not 
become significant until holes dominate and provide the conducting path. One can level the 
onset of conduction (or the apparent onset of conduction, as will be discussed later; ) as the 
threshold voltage V t h and define it as where the surface potential goes through intrinsic to 
a value of o s = 2 tpf. The threshold voltage is thus the gate voltage that produces a gate 
charge which just counterbalance the charge contained in surface states and in a depletion 
region that supports a voltage of 2 <pf. Further increases in gate voltage past V th result in an 
increase in the number of holes, thus enhancing conduction between the source and drain. 

It must be kept in mind that the start of actual conduction is not an abrupt process, in which 
channel is completely depleted of carriers at a given gate voltage Vq and then immediately 
acquires a finite inversion layer at a small increase in the magnitude of V G . The minority 
and majority carriers within the channel are increasing and decreasing continuously and at 
a finite rate for changes in the gate voltage. Goldberg described the situation very nicely 
when he said "The transition from depletion to inversion is a continuous process, and one 
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must be aware of the fact that the minority carrier concentration increases as the majority 
carrier concentration decreases and that carriers are always present.” 

Examination of the gate capacitance as a function of the gate voltage can yield additional 
insight into the physical operation of the MOS. Starting with the accumulation condition, 
figure A. 2 shows that the metal gate and the accumulation layer at the silicon surface firm 
a simple parallel-plate capacitor with a separation distance of t oa; and a normalized value 
of C /C 0 = 1. Capacitance is measured at the gate, with the source, drain, and substrate 
grounded. As negative bias is applied, the depletion region that forms tends to separate the 
two plates of the capacitor and reduces the capacitance. (The two plates are now the gate 
electrode and the bulk silicon at the edge of the depletion region. Separation is t ox + xq). 
Further application of a negative gate bias attracts holes to the surface (inversion), which in 
effect decreases the distance seperating the two plates, thus increasing capacitance. 

A. 3 Quantitative Analysis 

1. Linear Region, figure A. 3 shows an idealized device together with the coordinate system 
that is used in the analysis. A brief outline of the ’’plan of attack” used to derive the device 
equations is as follows: 

1. The channel-current density is integrated over a cross section of the channel (W dx) to 
obtain the current. 

2. Channel current now a function of the charge in the channel. 

3. This charge is found by summing all the system charge to zero. 

4. Charge is now related to the gate voltage by the use of Gauss’law. 

5. The expression for channel current can now be integrated over the length L of the channel. 

6. Channel current can now be equated to the external terminal current. 

Channel current can be written as 


I C = WJ J c (x,y)dx (A.l) 

Where W indicates channel width in the z direction. W is perpendicular to the direction of 
current flow. 

From Ohm's law, 
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P-chorinel MOS device 





Diagram* illustrating the position of charge under the gate as a function of bias: (a) accumula- 
tion case; ( b ) depletion case; (c) inversion case; ( d ) C vs V q plot of MOS gate showing the three 
characteristic regions together with a plot of x/\Td\ vs. V q. 


Figure A. 2: 





(b) 

(a) Diagrammatic .ketch of an MOS showing dlmen.lon. and direction, used for analy.li. 
(6) Diagram lllu.tratlng the Intrlnile and extrinsic element, of an MOS. 


Figure A. 3: 
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.A ( x j y'j uijx^Ey QfipJP (x) Ey 


Ic — W J q/j,pE y P(x)dx 


(A. 2) 


(A-3) 


Ic = Wq/j,pEy J p(x)dx (A. 4) 

i p is constant and independent of x . p p is a positive number, while p t - is negative.] 
d\v E y = -(^ 7)1 so that 


r r// dV 
-I c = Wpp—q 
dy 


J p(x)dx 


(A.5) 


fheveq f p(x)dx represents the mobile charge per unit area in the channel. The problem 
now reduced to evaluating q f p(x)dx. Since the total MOS-system charge must equal zero 


Qg + Qss + Qc + Qd — 0 (A. 6 ) 

/hero 

'c+Qss represents all the charge outside of the semiconductor material proper and Qc+Qu 
'presents all the charge within the semiconductor material . The channel charge is thus 


Qc — Qg + Qss + Qd 


(A.7) 
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The charge induced by the gate can be related to the gate voltage by Gauss’s law , which 
states 


f E ox (S)dS = 

J tax 


(A.8) 


dy — Q total 


(A.9) 


E 0 x is defined as 

ay,* 

dx 

where d\ ' ox = voltage across the oxide 
dx= oxide thickness now 

d Vox _ _ A Vox 

dx Ax 

Where 

A.x — "l - t ox 

AV 0X = -\V G - V(y)] 


(A.10) 

(A.ll) 

(A.12) 


so that 


dV ox V G -V(y) 
- + 


dx ' t ox (A ‘ 13) 

(The voltage across the oxide is simply the gate voltage minus the voltage on the channel, 
channel voltage will be a function of the distance in the y direction, ranging from V D at the 
drain to zero volts at the source.) Inserting Equation A. 7 into Equation A.13 and letting 


tox 


= C( capacitance per unit area ) yield 


Qc = [V G - V(y)]C (A. 14) 

Equation A. 14 relates the product of the gate capacitance per unit area and the voltage 
across the oxide to the charge per unit area under the gate. The gate charge found : 


Qc = _[y G _ V (y)]C - (Qss + Qd) (A.15) 

Equation (A.13) is a mathematical statement of the amount of mobile charge contained in 
the channel per unit area. Keep in mind that this is the charge that allows conduction 
between the source and drain. Qc will be enhanced by the gate voltage V G , decreased by 
the channel voltage created by the drain supply, V(y) , decreased by the charge stored in 
the depletion region beneath the channel, Qo, and either increased or decreased by Qss , 
depending upon its polarity. Later it will be interesting to investigate the condition required 
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the cause the channel charge to go to zero. 

Note that Qc — Q f p(x)dx, which is the desired quantity of Eq.(A.5). Substituting Eq.(A.5) 
into Eq.(A.13) yields 


Factoring out C, 


-Ic = Wy p --[V G - V(y)]C - ( Qss + Qd) 
ay 


(A.16) 


I c dy = Wy p CdV{ V G - V(y)) + 


Qss + Qd \ 


C 


7 


(A.17) 


Equation A.17 can now be integrated between 0 and L for the length and between 0 and Vp 
for the voltage: 


Icjdy = Wfx p C[V G £ D dV- I V(y)dy + 9 h + Qd J dy] (A.18) 

I C L = Wy p C(V G V D - 1/2 Vl + 9 . s ? . ± . 9p . V d ) (A.19) 

Now 

C = C 0 /WL (A.20) 

and 

Ic = -~^[~V G - V th )V D + 1/2V*} (A.21) 

where 

V t h = — ( Qss + Qd)/C. (A.22) 

(the threshold-voltage term will be discussed in more detail in 
Equation A.21 can also be written as 

Ic = -e;-(V G - Va.il'n + 1/21'=] (A.23) 

0 = (A.24) 

AjLqx 

Where (Because of the sign convention for p,/3 p is positive and /?„ is negative.) Now the 
channel current of the Equation Ic + Id — 0 orI G = —Id- The final form of the equation for 
the drain current of a device in the triode region can be written as 


(A.25) 


Id = -Wye - V a )V D - 1/2FJ] 
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This equation A. 25 assumes no drain or source resistance . At this point, it is- easy to in- 
clude the effects of a given parasitic resistance Rd and R s . Equation A. 17 is integrated from 
0 to L on the left hand side and from to on the right hand side. This integration yields 

Id = P[(V G - V th ){V' D - V' s ) - 1/2(V£ 2 - Vf )] (A.26) 

Which represents the triode region drain current. 


Vi, = V D -I D R D 


V's = I s Rs 


The saturation drain current can be written as 


Id = ~(3 


(Vg Z v p ) 2 

1 - 0R d (V a - V p ) + /il^20RAy^A^)) 


(A.27) 

(A.28) 


(A.29) 


Equation A.29 was used to generate the triode portion of the charactristic curves for P- 
channel device . 

The output current in the triode region depends upon the drain voltage, and hence the 
drain source terminals are like a resistor whose value depends upon the gate bias. Using 
a small signal approximation, the resistance is linear and adjustable. Close to the origin 
and when \Vp\ « [ 1'g ~ V t h | , the resistive characteristics are quite linear over a range of 
voltages and currents. Since the curves pass through the origin, the characteristics suggest 
that the MOS can be used as a voltage controlled d-c, as well as a-c, resistor. Equation 
was derived on the basis of the gradual channel approximation and is valid only for the case 
where 1 1 Id ( << He — Ut/>|- When \V D \ >= \V G — V p \, the device enters the saturation or 
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pinch off region . Equation A.29 is not also valid in the lov. r current region when Vq = V t h- 
Here the inversion layer is not yet fully formed, and the gate electric field terminates on 
comparable quantities of mobile charge (in the channel area) and immobile charge (in the 
depletion region beneath the channel). Because of this fact, the gate voltage no longer 
controls the channel conductance in the manner stated above equation. The lower boundary 
on the gate voltage such that this equation is still valid occurs where the channel carrier 
concentration is just equal to the bulk doping; that is, p C hannei = Np>. This occurs when 
the intrinsic Fermi level has been bent such that 4> s = 2 4> f . This subject will be covered in 
more detail in Saturation region. For a given voltage such that |V G | > \V th \, as the drain 
voltage is increased in magnitude from zero volts, the drain current increases linearly first, 
then slows down, and finally tends to level out as \V D \ is made large. The leveling out of 
the drain current is associated with the pinching off of the channel near the drain. Once the 
channel has pinched off, the current is said to have saturated at a given level and is then, 
to a first, approximation , independent of drain voltage. Pinchoff occurs because the voltage 
across the oxide falls below a critical value. The channel IR drop is the factor causing the 
reduction in electric field. When E ox is decreased to such a that it can not support sufficient 
mobile charge in a given portion of the channel, then that region decreases to approximately 
zero thickness and is said to have pinched off. There are two methods for approximating the 
mathmatical boundary between the triode and saturation regions. The first method involves 
letting the charge in the channel go to zero. This can be seen mathematically. Because the 
channel voltage 

[V G -V(y)]C = Qss + QD (A. 30) 

is highest at the drain, pinchoff will begin at the drain as |Vfi| is increased. V(y) in equation 
A.30 can thus be replaced by V D . Rearranging the terms of equqation gives the relationship 

V D = V G - V P (A.31) 

where 

Vp — —(Qss + Qd)/C. (A. 32) 

Notice that pinchoff voltage given is exactly the same as the threshold voltages . Equation 
A.32 states mathematically the boundary line between the triode and saturation regions 
which is shown plotted in figure. To the right of this line, operation is in the triode region, 
where |lbl < | Vg ~ 
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The second method of defining the saturation triode boundary can be found by examining 
equation for Id- For a given gate and threshold voltage, the magnitude of the drain current 
increases as |I.d| is increased from zero volts. Initially, the (Vg — Vh) Vd term dominates 
the expression, resulting in the increase of \Id\- Soon, however, the quadratic term becomes 
significant, with the result that the rate of increase of the current falls off. At some drain 
voltage. \Id\ will reach a maximum. Past this point, the equation predicts a decrease. How- 
ever, at the maximum current point, the device has reached saturation and the device model 
and equations change. This is why a decrease in current is not seen. To find the maximum 
of \Id\- the expression can be differentiated with respect to the drain voltage and set equal 
to zero. 

Once saturation has been reached, the voltage drop across the inverted portion of the channel 
tends to remain fixed at Vq — Vp, while Vp varies. To a first approximation, this constant 
voltage across a constant channel resistance results in a constant drain current. Once satura- 
tion has been reached . the output characteristics curves can be approximated by horizontal 
lines. The equation for the current in saturation can be found by placing V D = V G — Vp into 
previous equation for I d , which results in 


I D = •“( Vq - Vp ) 2 

This equation is valid for 

|Vd| >= \V G - V P \. 


(A. 33) 
(A. 34) 


A. 4 Ion- implanted Depletion mode devices 

Depletion mode MOSFET's are fabricated by introducing impurity ions of the opposite type 
to the starting substrate to form a shallow layer underneath gate of the device. This layer 
is most commonly achieved by ion implantation because excellent control of the number of 
impurity atoms and their distribution is possible to produce an n-channel depletion mode 
device in a standard NMOS process. N-type of impurities are selectively introduced into 
the P-type starting material. If an n-doped polysilicon gate is used, the resulting device 
will have a negative threshold voltage and is capable of conducting current with zero gate 
to source voltage. For this reason, these devices are typically used in logic circuits where a 
constant current load is desired. ' - 
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B.l Effects Of The DOS On The I-V Characteristics 
Of POLY-Si TFT / 

A high DOS results in a slow turn-on of the transfer characteristics. This makes difficult 
to define a threshold voltage (V T ) and a field effect mobility (/u), since a linear I D — V GS 
regime may not be present at all. The apparent threshold voltage (V' T ) and mobility (p) are 
deduced here from the Id=0 intercept and slope, respectively, of the tangent to the G — Vqs 
curve. The effect, of the DOS on the G — V GS characteristics is studied in detail by using the 


following set of equations, 

dG n(il> s ) — n Q 

dips ~ F s 

(B.l) 

F 2 = f Ec N(E)H(ip s ,E,E F )dE 

€ J E/ V 

+.4 f°° (E - E c ) 1/2 H(ip s ,E, E F )dE 

J S c 

(B.2) 

1 gs — f fb = tps + V ox = tp s + —(eF s — Q ss ) 

£ CX 

(B-3) 

N(E). in the upper half of the forbidden gap 


N(E) - At exp ( ^ ) + N D ex p( ^ ) 

(B.4) ; 

Expression for G, 

G = - Vfb - (A + 

(B.5) : 

The threshold voltage is defined as : 


Vr — tp s + t ox + Vfb 

(B.6) 


67 
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associated to the grain boundary defects are considered as uniformly distributed over the 
grain \olume. The band bending in the polysilicon is determined by Poisson’s equation: 


cPxj; q 

-fat = -Wtr°p + n(x) - no] 


(B.XO) 


Where N trap is the net concentration of negatively charged centers, and n(x)-no is the change 
in electron concentration with respect to equilibrium. These quantities can be written as 
functions of ib: 

Ntra P = f ‘ N(E ) [/(£, E f + #) - f(E, Ep)]dE 


n(ip) -n 0 = A [°°(E- E c f' 2 {f - f 0 )dE 

J JBc 


where: 


2N C 


f{E , E f + qip) = [1 + exp( 


7 r i/2fcT 3 /2 

E — E c — qip 


(B.ll) 
(B.12) 

(BJ3) 
(B-14) 

fo = f(if> = 0) (B.15) 

Multiplying both sides of Equation B.2 by 2^ and integrating from x=0 (oxide/semiconductor 
interface) to x=d (the unmodulated neutral bulk where ^ = 0)one obtain: 

2q 


kT 


)]-\ 


and 


F s 2 = ^ r dip [ ° N(E)[f(E, E f + qib) - f(E , E F )\dE 

€ JO JE V 


where 


+—A f V " dip H (E - E c ) 1/2 [f(E , E f + qib) - f(E , E F )dE] 

6 JO J JBc 

dip 


F, ~ dx 


(B-16) 


(B.17) 


at x=0 is the surface electric field. By changing the order the integration, This equation 
B.16 is reduces to: 


F 2 = ^ [ Ec N(E)H(iP s ,E,E F )dE 

6 JE V 

+— [°° (E-E c yt 2 H(ip z ,E,E F )dE 
€ J E c 

where 

E, E F ) = ^Zn(/ 0 [exp(^) + exj>( E j~r* )]) - ip SJ 
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The conductance G is defined as: 


G — qHn / n 

Jo 

By changing the variable of integration from x to ip and taking 
sides with respect to ib s ■ 



(B.20) 

the derivative of the both 


dG _ n(ip a — fto dG 
dib s q ^ n F s dijj s 

Gauss's Law can be written as: 


= QHon 0 


exp (qip s /kT) - 1 
_ 


(B.21) 


Vgs — Vps — ibs + V ox — ip s + —{eF s + Q ss ) (B.22) 

€ox 

Here, Q ss is the charge in the interface states at the oxide/semiconductor interface. This 
Equation (B.22)is re-written in terms of the fixed charge {Q fixed)-, by using the relationship: 


1 r d 

F s — - / (qn(x) + p sc )dx 

6 JO 

F s = f - f p sc (x)dx 

ep. n t Jo 

= G_ + Qsc 

ep n e 


(B.23) 


B.3 Determination of the DOS by the temperature 
Method 

Taking Equation (B.21), for qu’/kT >> 1 , and multiplying both sides by dip s /V GS , one 
obtains: 

dG _ Go exp(qip s /kT) 
dV^~~T F,(dV G /dv a ) 

This indicates that v 5 (V’ G5 ) can be deduced from a plot of log (dG/dip s ) vs 1/T. 


(B.24) 
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C.l Program for MOSFET Model 

#include <stdio.li> 

#include <math.h> 

double I_d() ; 

double Il_d() ; 

double vds [25] , value; 

double a() ; 

double gamma () ; 

double gammal () ; 

double l_i(); 

double v_b_s=-3.0; 

double v_g_s=5.0; 

double v_t=.045; 

double theta () ; 

double thetalO; 

double etalO; 

double eta() ; 

double I_d_sat=l . Oe-3 ; 

double k() ; 

double W=2. 16e-6; 

double velocity_l=15 . 0e4; 

double v_d_sat(); 

double v_c(); 

double a() ; 

double A() ; 

double B(); 

double q= 1.6e-19; 

double n=1.2e21; 

double d=10.0e-9; 

double D_j=0.5e-6; 

double epsilon_o=8.850e-12; 

double epsilon_s_i=12.0; 
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double 1=13. 4e-G; 
int i; 

double v_t2=4.0e-2; 
double v_tl=8.0e-2; 
double ph.i_f=0.3; 
double v_b_s2=-8 . 0 ; 
double v_b_sl=-l . 5 ; 

double ga.mma_0=49 . 5e-6 , t_o_x=0 . 045e-6 ; 


double I_d(i) 
int i ; 

{ 

double I_d; 

/♦printf ("@@@@7,lf\n" ,vds [i] ) ; ♦/ 

if ( (vds [i] >=0 && vds[i]<v_d_sat()) && v_g_s >=v_t) 

{ 

I_d= gamma(i) *W* (v_g_s-v_t-a() *vds [i] /2 . 0) *vds [i] /I ; 

/*printf (" ! \t*/.lf \n" ,gamma(i) ) ; 
printf ("5\t7.1f\n" , W) ; 
printf ("t\t7.1f\n" ,v_g_s) ; 
printf ("r\t7.1f\n" ,v_d_s) ; 
printf C‘r\t7.1f\n" ,v_t) ; 
printf ("r\t7.1f\n" ,a() ) ; 
printf ("r\t7.1f\n" ,v_d_s) ;♦/ 

/♦printf ("I_d=7.1f \n" ,I_d) ;*/ ■ 

} 

else if (vds [i]>v_d_sat () && v_g_s>=v_t) 

{ 

I_d=gamma(i) *W* (v_g_s — v_t-a() *v_d_sat () /2 . 0) *v_d_sat C) /l_i (i) , 
/♦printf ( “ganma ! \t7.1f \n" , gamma 0 ) ; 
printf ("WYtyilf \n" ,W) ; 



appendix c. 


74 


printf ("v_g_s\t*/.lf \n" , v_g_s) ; 
printf ("v_d_s\t7.1f\n" ,v_d_s) ; 
printf ‘/.If \n" ,v_t) ; 

printf ("a\t*/.lf\n" , a() ) ; 
printf ("v_d_sat()\t7.1f\n" ,v_d_satQ) ; 
printf C"l_i()\t'/.lf\n" ,l_i()) ;*/ 
/♦printf ("I_d=7.1f \n" , I_d) ; */ 

> 

else if (v_g_s<v_t) 

{ 

l_d=0 . 0 ; 

/♦printf ("I_d=7.1f\n" , I_d) ; ♦/ 

Jreturn I_d; 

> 


double mu_o(gamma_0,t_o_x) 
double gamma_0,t_o_x; 

{ 

double epsilon_o=8 . 85e-12; 
double epsilon_o_x=3 . 9 ; 

return (gamma_0*t_o_x) / (epsilon_o*epsilon_o_jJ , 
return ; 

> 


double v_c() 

return (velocity_l*l/mu_o(49.5e — 6,0.045e 6)), 
/♦printf ( " v_c\t7.1f \n" , v_c 0 ) ; 

printf ( "mu_o\t7.1f \n" ,mu_o (49.5,0. 045e-6) ) ; 
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printf ("velocity_l\f/.lf \n" , velocity.].) ; 
printf ("lengthYt/ilf\n" ,1) ; 

*/ 

> 

double k() 

{ 

return (v_u2-v_t 1) / (sqrt (2 . 0*phi_f ■ -v_b_s2) -sqrt (2 . 0*pbi_f-v_b_sl) ) ; 

> 

double v.d.satC) 

{ 

double V.ds , suml , sum2 , sum3 ; 
suml=x r _g_s-v_t ; 
suml=suml/l . 0; 
sum2=v_c() ; 

sum3=suml*suml+v_c()*v_c() ; 
sum3=sqrt (sxrm3) ; 

V_ds= ( ( (v_g_s-v_t ) /a() ) +v_c () ) -sqrt ( ( (v_g_s-v_t) * (v_g_s-v_t) /pow(a() ,2.0)) 
+(pow(v_c() ,2.0))) ; 

V_ds=suml+si:m2-sum3 ; 

/ *printf ( " V.Dst aurated=\t'/.lf \n" , V.ds) ; */ 

return V.ds ; 

> 


double a() 


double temp=0.0; 
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temp=l . 0 ; 

temp+=0.5*kO*(1.0 - 1.0/Cl. 41 +0.43*(2.0*phi_f - v_b_s)))/ 

sqrt(2.0*ph.i_f-v_b_s) ; 
return temp; 

} 


double gamma (i) 
int i; 

{ 

double gamma_0=49 . 5e-6 ; 
double t ; 

/*printf ("##7.f\n" ,vds[i]) ;*/ 

t = (gamma_0/ ( 1 . 0+theta ( ) * (v_g_s-v_t-vds [i] /2 . 0+k () *sqrt (2 . 0*phi_f-v_b_s) )+ 

eta()*vds [i] )) ; 

/*printf ("*/.lf \n" ,t) ; */ 
return t ; 

} 


double thetaO 

{ 

double theta_0=0 .055 ; 
double theta_l=0 .004e-6; 

/*printf{ "theta* y.lf\n" ,theta_0+theta_l/l) ;*/ 
return theta_0+theta_l/l ; 

> 
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double eta() 

{ 

double eta_0=0.008; 
double eta_l=0.25e-6; 

/*printf ("eta=%lf\n" , eta_0+eta_l/l) ;*/ 
return eta_0+eua_l/l ; 

> 


double l_i (i) 
mt i; 

{ 

return 1-1* (sqrt (pow(v_c() ,2 . 0)+2 . 0*A()* (1 . 0+B()*I_d_sat)* (vds [i] -v_d_sat () ) ) 
/(A()*(1.0+B()*I_d_sat)); 

} 

double A() 

{ 

return q*n*l*l/ (epsilon_o*epsilon_s_i) ; 

> 

double B() 

{ 

return 2.0*(log(D_j/d)-l .0)/(q*n*velocity_l*W*D_j) ; 

> 
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double gammaK) 

{ 

double gajnmal_0=49.5e-6; 

return g amm a 1 0/ ( 1 . 0+thetal () * (v_g_s-v_t-vds [i] /2 . 0+k() *sqrt (2 . 0*phi_f-v_b_s) ) +etal () 
> 

double thetalO 

{ 

double theta_0=0 . 055; 
double theta, 1=0. 004e-6; 
return theta_0+theta_l/l_i (i) ; 

> 

double etalO 

{ 

double eta_0=0.008; 
double eta_l=0.25e-6; 
return eta_0+eta_l/l_i (i) ; 

> 


mainO 

{ 

int i ; 

FILE *data; 

data=f open ("data. dat_ .055_5" , "w") ; 

printf ("the saturated, voltage=\t*/.lf\n" ,v_d_sat()) ; 
printf ("enter 25 real numbers\n") ; 
for(i=0; i<25;i++) 
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-r 

scanf ("7.1f " ,&value) ; 
vds [i] =value ; 

} 

/*printf ("PLEASE ENTER v_t , v_g_s , v_d_s , threshold voltage,voltage_drain_sourc 

voltage_gate_source\n") ; 

scanf ("’/.If '/.IfXlf" ,&v_t ,&v_g_s ,&vds) ;*/ 

/* printf ("cite gamma=\t*/.lf \n" , gamma (i) ) ; 
printf (“the v_c=\t%lf \n" , v_c()) ; */ 


f or(i=0;i<25 ; i++) 

{ 

fprintf (data, "7.1f \t*/,lf\n" , vds [i] ,I_d(i)) ; 
printf ("vds=%Lf\n" , i , vds [i] ) ; 
printf ("I_d[‘/.cj=7.1f\n",i,I_d(i)) ; 

/* I.dCi);*/ 

> 

printf ("\n") ; 
f close (data) ; 

> 



APPENDIX C. 


80 


C.2 Program For The Subthreshold Characteristics Of 
Poly-Si TFT’s 

tinclude <ma-th.li> 

#include <stdio.h> 
double y() ; 
double x; 
double g_0(); 
double epsilon_ox=3 . 9 ; 
double t_ox=0 . le-6; 
double n_t(); 
double n_c=l . 3e25 ; 

/♦double n_c=3.2e25;*/ 
double (e_fec)=-0 .42; 
double mu=5e-4; 
double q=1.6e-19; 
double t; 
double v_g; 

double k=(69 . Oe-3) /800 . 0 ; 
double k_0(); 
double G() ; 
double t_g=800 . 0 ; 
double v_0(); 
double k_tg=(69.0e-3) ; 
double n_g=3.0e26; 
double d=0 . 6e-6 ; 
double f 0 ; 

double epsilon_0=8.85e-12; 
double epsilon_s=ll .9; 
double y() 

return loglO (g_0 C) *epsilon_ox*epsilon_0/ Cn_t 0 *d*t_ox) ) + (2.0*t_g/t- 2 . 0) * (logl 
-logl0(v_0())) ; 
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> 

double f() 

{ 

return sqrt (2 . 0*n_t () *k_tg*q/ (epsilon_s*epsilon_0) ) ; 

} 

double n_t() 

■C 

return n_g*k_tg*exp ( (e_f ec) / (k_tg) ) ; 

> 

double g_0() 

{ 

return mu*d*n_c*exp((e_fec)/(k*t)) ; 

> 

double v_0() 

{ 

return epsilon_s*t_ox*f () /epsilon_ox ; 

> 

double GO 

return g_0 O+k_0O *pow(x,2 . 0*t_g/t-l .0) ; 

> 

double k_0() 

{ 

return q*mu*n_c* (t / (2 . 0*t_g-t) ) *pow(2 . 0*epsilon_s*epsilon_0/ (n_g*pow(q, 2 . 0) ) ,0 . 
pow (pow (epsilon_ox , 2 . 0) / (2 . 0*eps ilon_s*epsilon_0*pow (t_ox ,2.0)* (n_g/q) *pow (k*q , 
pow(t_g,2.0)) ,2.0*t_g/t-l .0/2.0) ; 

> 

mainO 

■C 

int i; 

FILE *data; 

data=f open("a.dtft4" ,"w") ; 
for (i=0;i<12;i++) 
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{ 

/♦printf ("Please Enter x,t\n") ; 
scant ('"/.If %lf" 

printf ("Please Enter v_g,t\n") ; 

scanf ( '"/.lf/.lf " . &v_g , &t ) ; 

printf ("y=\ty.lf\n" ,y (x)) ; 

printf ("x=\ty.lf\n" ,x) ; 

printf ("x_2=\ty.lf\n" , x*x) ; 

printf ("t=\ty.lf \n" ,t) ; 

printf ("f=\t%lf \n" ,f ()) ; 

printf ("n=\tyif \n" ,n_t () ) ; 

printf ("g=\ty.lf\n" ,g_0()) ; 

printf ("v=\ty.lf\n" ,v_0O) ; 

printf ("k=\ty.lf\n" ,k_tg) ; 

printf ("n_g=\t*/.lf\n" ,n_g) ; 

printf ("n_c=\t*/.lf \n" ,n_c) ; 

printf ("e_f=\t'/.lf\n" ,e_fec) ; 

printf ("*/.f\ty.lf\n" , (x*x) ,pow(10 ,y() ) ) ; 

/♦fprintf (dana, '"/.t \t*/.lf \n" ,loglO(x*x) ,y()) ;*/ 

/*f printf (dana, "*/.f \t'/.lf \n" ,loglO(x) ,y()) ;*/ 

fprintf (data, "*/.f\ty.lf\n" ,log(G() ) ,v_g) ; 

> 

f close (data) ; 

} 
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C.3 Program For Accumulation Mode Of Poly-Si TFT’ 

#include <stdio.h> 

#include<math . h> 

double mu_b=ll . Oe-4; 

double q=1.6e-19; 

double mu() ; 

double psi(); 

double W=10.0e-6; 

double L=10.0e-6; 

double psi_0=0.55; 

double alpha () ; 

double I () ; 

double V_d=0.5; 

double V_g; 

double N_st=2.1el6; 

double epsilon=8 .85e-12; 

double C_ox=12 . 0* (8 . 85e-12) / (78 . Oe-9) ; 

/★double C_ox=12 . 0* (8 . 85e-12) / (78 . Oe-9) = . 0013615385 ; */ 

double d=10.0e-9; 

double V_th=l .2; 

double k=(8.614e-5)*(1.6e-19); 

/★double k= (8 . 614e-5) * (1 . 6e-19) =1 . 37824e-23 ; ★/ 

double T=300.0; 

void input .variables () 

{ 

printf ("ENTER V_g\n") ; 
scanf ("'/.If " ,&V_g) ; 

> 

double mu() 

{ 

return mu.b+exp (-q*psi () / (k*T) ) ; 

> 

double psi() 
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{ 

return alphaO/ (V_g-V_th+alpha()/psi_0) ; 

> 

double alpha O 

{ 

return q*q*N_st*N_st*d/ (8 . 0*epsilon*C_ox) ; 

> 

double 10 

return (W/L)*C_ox*mu()*( (V_g-V_th)*V_d-pow(V_d,2.0)/2.0) ; 

> 

main() 

{ 

int i; 

FILE *data; 

data=fopen("data. datanuli3" , "w") ; 
f or (i=0 ; i<25 ; i++) 

■C 

mput_variables () ; 

printf ("p=\t*/.e\n" ,psi()) ; 

printf ("v=\t7.If\n" ,V_g) ; 

printf ("c=\t7.Lf\n" ,C_ox) ; 

printf ("k=\t*/.e\n" ,k) ; 

printf ("kT=\t*e\n" ,k*T) ; 

printf ("ex=\t*>.e\n" ,exp(-q*psi()/ (k*T))) ; 

printf ("q=\t*/.e\n" ,q*q) ; 

printf ("I=\t'/.e\n" ,10) ; 

printf ("ql=\tXe\n" ,q*N_st) ; 

printf ( " V_=\tXe\n" , pow (V_d ,2.0)); 

printf C"a=\t*/,e\n" , alphaO) ; 

printf ("n^N^/.eXn" ,mu() ) ; 

printf ("ml\t7,lf\n" ,mu()) ; 

f printf (data, "Xe\t7.e\n" , V_g , 1 0 ) ; 
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f close (data) ; 

> 
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